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We have studied the critical temperature of Diluted Magnetic Semiconductors by means of Monte 
Carlo simulations and Coherent-Potential- Approximation (CPA) calculations. In our model for 
this system, the magnetic ions couple with the carriers through an antiferromagnetic exchange 
interaction, J, and an electrostatic interaction W . The effective impurity potential J — W controls 
the hybridization between the magnetic impurities and the hole charge on the dopants. We find 
that the critical temperature depends substantially on the hole charge on the magnetic impurities. 
The CPA critical temperature is always lower than the obtained in the Monte Carlo simulations, 
although all trends in the simulation results are reproduced in the CPA calculations. Finally we 
predict the existence of pockets of phase segregation instability close to the carriers band edges. 

PACS numbers: 75.50.Pp, 75.10Lp 
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I. INTRODUCTION. 

It has been recently discovered that some Mn 
doped semiconductors, such as Gai-^Mn^As, undergo 
a ferromagnetic-paramagnetic transition at temperatures 
near 100K [Q, gj. Motivated by this high Curie temper- 
ature, there have been further experimental studies in 
III-V magnetic semiconductors, and ferro-paramagnetic 
transition temperatures as high as 1000K have been ob- 
tained in (Ga,Mn)N gj. The optimal Mn concentration 
is near 5% in these materials, making magnetic ions to 
be rather diluted in the host semiconductor. 

The existence of ferromagnetic semiconductors at room 
temperature opens the possibility of combining the spa- 
tial modulations of both carrier and spin density. There- 
fore, diluted magnetic semiconductors (DMS) have at- 
tracted much attention for their potential use in spin- 
tronic devices. Some spintronic devices have been al- 
ready implemented [Q, and different heterostructures in 
order to optimize and control the critical temperature 
have been proposed ^ |[. 

In this work we study magnetic properties III-V mag- 
netic semiconductors, in particular of GaMnAs. The high 
critical temperature DMS's have a high concentration of 
randomly located manganese ions. From electron para- 
magnetic resonance and optical experiments |)| [To), it is 
known that the Mn ions in these semiconductors have 
a S = 5/2 local moment. The Mn ions go in substi- 
tutionally at the cation sites, Ga, and contribute holes 
to the semiconductor valence band. Although the lo- 
cal moment of the magnetic ions in GaMnAs indicates 
that each Mn should contribute with a hole, it is exper- 
imentally observed that the density of itinerant holes, 
p, is a small fraction (~ 10%) of the Mn concentration 
jDJ . The reduction in the carrier density is probably due 



to the presence of antisite defects which, acting as deep 
donors partially compensating the Mn acceptors, reduce 
the hole density jH| . A large number of antisite defects 
is expected to be present in these systems because, in 
order to minimize the mobility of the Mn ions and the 
existence of clusters, the semiconductors are fabricated 
using molecular beam epitaxy growing techniques at low 
temperatures, for which a large number of these defects 
is known to occur 12 T|. 



It is widely accepted that the ferromagnetism in III-V 
semiconductors is induced by the motion of the electrical 
carriers in the system. However, there is no agreement on 
the precise role of the carriers. They can be either itin- 
erant, with a clear character of the host semiconductor 
valence band states jl4|, |l5| [if], |l7j , or rather they can be 
mostly located on Mn impurities, with a clear impurity 
band magnetic, polaron character |l8|, [H], ^lj. On 
the other hand, the carrier induced ferromagnetic inter- 
action between Mn 2+ ions in magnetic semiconductors 
has some resemblance to the one observed in manganese 
perovskites p2| , double perovskites Q , and pyrochlores 
P4j . In these systems, the phase diagram is different 
from that of conventional itinerant ferromagnets, show- 
ing first-order transitions, phase separation |2q| , and for- 
mation of magnetic polarons near the Curie temperature 

In order to study these issues, we use a tight-binding 
model Hamiltonian in which, besides the Hund's cou- 
pling, a site energy difference between magnetic and nor- 
mal ions, W, is introduced. Monte Carlo (MC) simula- 
tions are then performed for this model, with particular 
emphasis on the W dependence of the phase diagram. 
A mean field calculation within the coherent-potential- 
approximation (CPA) is then carried out and compared 
with MC results. The CPA treatment confirms the ro- 
bustness of the MC results for the effect of W, in addition 
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to allowing a simpler exploration of important issues such 
as the existence of phase separation. 
Our main results are the following: 

1) The hole charge on the Mn ions, controlled by the 
impurity potential W , has a profound effect on the tran- 
sition temperature. Upon increasing W, the hole charge 
on Mn decreases, leading to a substantial decrease of the 
Curie temperature. 

2) For much of the model's parameter space, the MC crit- 
ical temperature is higher than that obtained using the 
simplest Weiss mean field theory. The MC critical tem- 
perature is always higher than the corresponding CPA 
temperature. Nevertheless, all trends in the MC results, 
particularly the effect of the impurity potential, are re- 
produced in the CPA calculation. We believe this to be 
a robust feature. 

3) The CPA calculation shows the existence of small 
pockets of phase segregation instability close to band 
edges. 

The paper is organized in the following way. In sec- 
tion II we introduce and explain the model. Section III 
and IV describe the MC and CPA methods, respectively. 
In Section V we present our results, both MC and CPA, 
making contact with previous calculations and/or exper- 
iments. Section VI summarizes our work. 



II. THE MODEL. 



localization effects. With all this, the Hamiltonian for 
the holes is written as 



It is generally accepted JTJ, H, H, H, @ that 
the DMS's are governed by the following Hamiltonian: 



H = H h 



Hi, 



H, 



(1) 



Hh and Hi ons are the parts of the Hamiltonian describ- 
ing the holes and the Mn ions, respectively. The term 
Hh-ions represents the interaction between the magnetic 
ions and the carriers. 

Hh is the sum of the kinetic energy of the carriers, 
the hole-hole interaction energy, and the interaction of 
the holes with the disorder potential from randomness in 
the host semiconductor, in particular from the presence 
of antisite defects. In the actual magnetic semiconduc- 
tors, the carrier density is of the order of 10 20 cm~ 3 , 
and, for these values, it is justified to neglect the effect 
of the carrier-carrier interaction. In this work we do no 
take into account the interaction of carriers with antisite 
disorder. This contribution is probably important in ac- 
tual samples, but here we are mainly interested in the 
effect of the on-site energy difference and thermal fluc- 
tuations on the critical temperature. We describe the 
motion of the holes using a single one-band hole disper- 
sion, instead of the more appropriated six bands envelope 
function formalism ^8|. This one-band model is proba- 
bly insufficient to give a quantitative prediction for the 
Curie temperature, but it includes the main features of 
the hole system we want to study: the effect of thermal 
fluctuations, Mn-GaAs orbital hybridization, and charge 



H, 



3," 



(2) 
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where the sum runs over all first neighbors pairs on a 
simple cubic lattice, t is the tunneling amplitude, and 
Cf a creates a hole with spin a at site i. This model 
has a hole bandwidth of 12£, and the parameter t has 
been chosen to give a bandwidth of the order of that of 
GaAs (~T0eV), implying t ~0.84eF. The band edges 
of this model correspond to a parabolic dispersion with 
an effective mass m* = ti 2 /2a 2 t, being a the lattice pa- 
rameter of the simple cubic lattice. The value of a has 
been chosen to make a 3 the volume per GaAs unit, i.e. 
a=3.56A. With these parameters, we obtain an effective 
mass 77i*=0.36mo. 

The term Hi ons in Eq.(|l|) describes the direct antifer- 
romagnetic interactions between the magnetic moments 
of the Mn ions. Those are much smaller than the inter- 
actions with the carrier spins and, therefore, we neglect 
them. 

Finally, the term Hh-ions represents the coupling be- 
tween the electrical carriers and the magnetic impurities. 
There are two contributions to this term: i) An antifcr- 
romagnetic exchange interaction between the spins S/ of 
the Mn 2+ ions located at sites I and the spins of the 
itinerant carriers. This interaction produces a long range 
ferromagnetic interaction between the Mn ions. And ii) 
an interaction between the carriers charge and the po- 
tential arising from the Mn dopants: 



Hh-ions = -X ^ ^t,a a cc,pCl,l} + W ^ Cf^ C If , 



(3) 



Here the sum runs over the positions of the magnetic im- 
purities, and er arc the Pauli matrices. The origin of W 
is Coulombic, coming from the different electronegativity 
of Mn and GaAs atoms, and from the screening of the 
impurities by the electrical carriers. The quantity J — W 
governs the hybridization of the magnetic impurities with 
the host bands, and the excess of charge on the impu- 
rities with respect the host semiconductor atoms. For 
instance, in the case of a perfect ferromagnetic ground 
state, the electrical charge is uniformly distributed for 
J=W. The value of J in our model is connected with 
the exchange coupling J p d appearing in the continuum 
model H^] through J = SJ p d/2a 3 . The value of J p d is 
typically 0M-0.leVnm 3 @, §, |(| yielding J -1.4- 
2.8eV, which gives a value of J/t ~1.7-3.4. Local spin 
density calculations (LSDA) and supercell methods |3^| 
applied to GaMnAs predict a value of J rather larger 
than the experimental one, partially because the LSDA 
calculations overestimate exchange coupling strengths. 
Lacking reliable experimental information on the value 
of W, we consider it as a parameter with values between 
and J. Summarizing the model: we restrict ourselves 
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to the terms Hh and H^-ions, given by Eq. || and Eq. [|, 
respectively. 



III. DESCRIPTION OF THE MONTE CARLO 
ALGORITHM. 

We have performed classical MC simulations on the Mn 
core spin angles in order to calculate the temperature, T, 
dependent magnetic phase diagram. The simulations are 
done in a N x N x N cubic lattice with periodic boundary 
conditions, and Nmu=xN 3 randomly located magnetic 
ions. The number of carriers is taken to be Nh=pN 3 . 

The simulations are done as follows. First, we produce 
a realization of the disorder: the positions of the Mn ions. 
For this realization, we compute the expectation value of 
the hole spin orientation at every place, 



-<C 



At a given temperature, the Mn spins fluctuate thermi- 
cally on a fictitious magnetic field [[33| proportional to 
the carriers spin density, 



B/OO 



J 

s 



(5) 



We perform MC simulations on the Mn spin angles until 
the system is equilibrated, when we choose a particular 
core spin configuration, and recalculate the hole spin den- 
sity. This cycle is repeated until convergence in physical 
quantities is attained. Calculating the hole spin density 
requires the diagonalization of the Hamiltonian Eq. ([!]) at 
each step of the cycle, is evaluated by filling up the 
Nh lowest energy levels of the Hamiltonian. This can be 
done because the Fermi temperature is much higher than 
any other temperature in the system. For instance, for 
cc=10%, p=l%, and J=W =3i, the Fermi energy is near 
0.5i, while the Curie temperature is ~0.015t. 

The CPU diagonalization cost grows as N 6 , severely 
limiting the system's sizes under study. By diagonalizing 
the Hamiltonian in real space, it is not possible to analyze 
systems bigger than N—6 p5fl . In order to study 
bigger systems, it is better to use an adapted basis rather 
than the atomic localized basis. In the simulation, we use 
the basis obtained by diagonalizing the hole Hamiltonian 
for a perfect ferromagnetic alignment of the Mn ion spins. 
For the updating of the electron spin density at each 
cycle of the simulation, it is not necessary to include all 
the eigenvectors of the extended basis. We have found 
enough to work with a reduced basis formed by the N c 
lowest energy eigenstates of this basis for each carrier spin 
orientation. The value of N c depends on the particular 
values of the parameters. In general, we find that the size 
of the required basis increases with the value of J- W, i.e. 
it increases with the degree of localization of the hole 
charge on the Mn ions. Typically we can obtain good 
convergence with relatively small values of N c , most of 
the results presented in this paper have been obtained 



with N c =57. Note that, in the extended basis we use, 
the square of the wave functions is not spatially uniform, 
as it would happen in a Bloch or plane wave basis, but it 
is modulated by the spatial fluctuations of the potential 
felt by the carriers. Using this extended basis and an 
appropriate cutoff N c , we can perform MC simulations 
in systems with sizes up to N=16. 

In the simulations, we calculate the thermal averages 
of different quantities: the absolute value of the Mn ions 
spin polarization, 



M 



1 



SN Mn 



(0) 



the absolute value of the hole's spin polarization, 

i<i5>i>> (7) 



(4) and the average value of the hole charge on the Mn ions, 



777:<E< 



(8) 



being pi =< J^a Cf a d^ a > the electrical hole charge on 
site i. In the previous expressions, the index / runs over 
the Mn positions and i over all lattices sites. 



IV. CPA METHOD. 

We have performed mean field calculations for the 
present model within the framework of the Coherent- 
Potential- Approximation (CPA) JH, |7|, §§]. In re- 
cent times, this formalism has come to be known as 
the Dynamical-Mean-Field- Approximation (DMFA) [§9) , 



and we use both terms (CPA and DMFA) as synonymous 
expressions. Given the possible confusion from the use 
of different flavors and names under the common um- 
brella of mean field calculations, we devote this section 
to briefly describe our CPA formalism, pointing out the 
origin and nature of the approximations and, where ap- 
propriate, making contact with other approaches. 

Tracing out the fermions at constant chemical poten- 
tial and temperature, we end up with a probability dis- 
tribution for the core spins given by the following Boltz- 
mann weight: 



P({Si})««Bph9fi({Si})] 



(9) 



where fl ({S^}) is the carriers grand canonical potential 
for a given core-spin configuration. The exact distribu- 
tion P({Sj}) minimizes the following thermodynamical 
potential, considered as a functional over spin distribu- 
tions: 



F[V] = (n({S i })) + T(log(V)) 



(10) 



where averages are taken over core-spin distributions. 
This formulation offers a convenient starting point for 
approximate variational approaches. 
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The common basis of all mean field calculations is the 
minimization of Eq. [To] with a factorization ansatz for the 
spin distribution: 



V({S l }) = l[p(S l ) 



From the critical phenomena point of view, this fac- 
torization legitimizes the term mean field, implying that 
all criticality (if any) will be classical. Carrying out the 
minimization of Eq. [l^ with the only simplification of 
the factorization could be termed the absolute mean field 
approximation. Unfortunately, this procedure leaves us 
with the task of calculating the grand canonical poten- 
tial for fermions (i.e. the density of states) in a site- 
disordered background: a formidable problem requiring 
numerical treatments (see, for instance, Ref. |Q for such 
a treatment in the Double- Exchange model). 

To cope with the previous difficulty, additional ap- 
proximations are needed to obtain the average density of 
states for the disordered problem. We will describe two 
of them, namely, the Virtual-Crystal- Approximation and 
CPA treatment (see Ref. 38 for a review on disordered 
systems) . 

i) Virtual Crystal Approximation ( VCA) . The average 
density of states for the disordered system is replaced by 
that of the average Hamiltonian. That implies a trans- 
lationally invariant system with site energies given by 
the following average (a 2 x 2 matrix in spin space): 

£a,/3 = x(W + J^jf- ■ ffa.p)- The minimization of T is 
now straightforward, leading to the following well-known 
expression for the transition temperature: 



(12) 



where n a (\x) is the density of states per site and spin 
at the Fermi level (degenerate fermions assumed). This 
scheme is often referred to as the Weiss mean field ap- 
proach, and we will consider both expressions (Weiss 
mean field and VCA) as synonymous terms. 

ii) Coherent Potential Approximation (CPA). The 
CPA is a time-honored approach for the site-disordered 
"alloy" problem |3(| |37|. It is known to produce good 
results, particularly in three dimensions [p8[ . The aver- 
aged Green's function is obtained from a translationally- 
invariant, effective medium, characterized by a (site- 
diagonal and energy dependent) self-energy E a> ^(z) in 
the following manner: 



(G(z)) ~ G 



(z-h y 



(13) 



where (Dirac notation used to emphasize the one-body 
character of the problem) 



H CPA =-t^2 \h CT ) 0") u \ + \h a ) Z>a,p(z) (i, 0\ 



with sum over spin indices implied. Given a particular 
site-disorder distribution, the self-energy is obtained with 
a minimum scattering requirement: the substitution of 
an effective site by a real one should produce zero (local) 

That is: 



(11) scattering on the average 41 



(£,-£)[! 



o 



(15) 



where hats mean 2x2 matrices in spin space, and the 
random quantities £i are given by: 

_ J 0, 'probability = (1 — x) 

/3) - | w + jS±,o- a ^ probability = xp(Si) 

(16) 

As usual, the minimization of the thermodynamic po- 
tential can be cast in the form of a self-consistency equa- 
tion for the probability distribution: 



p (Si) oc exp [—(3 5f2(Si 



(17) 



where <5f2(Sj) can be interpreted as the total change in 
the fermion's grand canonical potential for embedding an 
impurity site in the effective medium: 



SQ(Si 



dufiu) — SlogdetM(Sj) (18) 

7T 



(14) 



f(tu) is the Fermi-Dirac factor, 9 means imaginary part, 
and Af (Si) is given by: 

M(Si) = 1 - gT(lo + ) [e(S.) - E( W +)] (19) 

It is through eqs. [l7| and [l8| that the CPA can be con- 
nected with the DMFA (see, for instance, Ref. Q), 
an approach primarily intended for genuine many-body 
problems like the Hubbard model j3jJ]. It has been 
recently applied to models of itinerant carriers with 
ion impurities, believed to be relevant for manganites 
and doped magnetic semiconductors (see, for instance, 
Ref. §§@ H). 

Before leaving this section, let us comment on some 
features of these approaches. It is clear that the VCA 
is simple and handy, but the approximation in the den- 
sity of states is too drastic. For instance, it predicts a 
ferromagnetic transition for any finite density of states. 
This cannot be true, even at the mean field level: for 
W = and small J, the problem becomes the standard 
Ruderman-Kittel-Kasuya-Yosida (RKKY) Hamiltonian, 
whose absolute mean field treatment shows ferromag- 
netism only for carrier density (holes or electrons) below 
a critical value p c = 0.5038 (sc lattice), due to the oscil- 
latory behavior of the spin-spin coupling. Furthermore, 
the VCA is not sensitive to the value of W, obviously a 
major limitation for this problem. In contrast, the CPA 
is free from these shortcomings. For instance, it is able to 
reproduce the absolute mean field results of the RKKY 
limit p6[ . In addition, extensive numerical results avail- 
able for the exact density of states in the closely related 
Double- Exchange problem [^tj , show virtually no differ- 
ences with the CPA density of states j|7| . This reinforces 
our confidence in this method for the present problem. 
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FIG. 1: Mn ions spin polarization as a function of temperature 
for J=3t, i=10%, p=0.7%, N=U and different values of the 
parameter W. Results for different realizations of disorder 
are plotted. 



FIG. 3: Average hole charge on Mn ions, as a function of T, 
for two different values of W . The results corresponds to the 
case J=3t, x=10%, p=0.7%, and N=14. The arrows mark 
the Curie temperature, as obtained from Fig. |l| and Fig. |^. 
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FIG. 2: Holes spin polarization as a function of temperature 
for J=3i, a;=10%, p=0.7%, JV=14 and different values of the 
parameter W . Results for different realizations of disorder 
are plotted. 



V. RESULTS. 
A. Monte Carlo Results. 

In Fig. we plot the Mn ion's spin polarization, M, 
as a function of T, for J=3i, and different values of W. 
The results correspond to a Mn concentration x=10%, a 
hole concentration p—0. 7%, and a unit cell size N=1A. In 
Fig. U we plot the hole's spin polarization for the same 
set of parameters. Due to the finite size effects (note 
that there are only 274 Mn ions) the polarizations M 
and m are different from zero at any temperature, and 
we define the critical temperature, T c , as the point where 
the second derivate of the polarization with respect to T 
changes sign. This way of obtaining T c implies uncer- 
tainties of around 5%. Within this error margin, the T c 's 
obtained from the temperature dependence of m and M 
coincide. This is the expected result as the carrier's spin 
polarization is due to the Mn polarization. 

Comparing Fig. |l| and Fig. ^, we see that the ion spin 



polarization decreases faster with temperature than that 
of the holes. This seems to be generally true for the range 
of hole and Mn densities of interest, and a value of J=3t. 
This behavior has been already found in Weiss mean field 
calculations || . Different disorder realizations are super- 
imposed in figures 1 and 2. We find that, for small values 
of J-W, different disorder realizations give practically the 
same polarization curves, whereas, for small values of W, 
fluctuations related with disorder become more impor- 
tant. This is because the larger J-W, the deeper the 
effective impurity level. This increases the charge local- 
ized in the Mn ions, making the fluctuations associated 
with impurity randomness more important. 

The critical temperature decreases when W increases. 
As we will show, this is related to the amount of hole 
charge localized on Mn ions. In Fig. ||we plot the average 
hole charge on Mn as a function of T, for two different 
values of W. Fig. || shows that the charge on the Mn's 
decreases as the effective impurity potential increases, 
that Qmu decreases with temperature, and that there is 
a change in the sign of the second derivate of Q Mn with 
respect to temperature at T c . 

The delocalization of charge with increasing T, and its 
abrupt change at T c , indicate a possible change in trans- 
port properties of the system near the ferro-paramagnetic 
transition: we expect that the electrical resistivity of the 
system will decrease when T increases. For a value of 
W=J, the effective impurity potential is exactly zero in 
the perfect ferromagnetic state. Close to this value, the 
charge on each Mn ion is practically p, independently 
of T . However, for values of W near zero, the charge 
on the Mn ions is very large, more than three times p. 
Note that, in the case of W—0, the critical value of J 
for binding a hole of the appropriate spin to a Mn ion is 
J c = 3.95677t. For values of J > J c , the holes form a 
parallel-spin impurity band split from the main general 
band ^5|. In that case, the physics becomes more similar 
to that of the Double-Exchange materials. In our model, 
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FIG. 4: Critical temperature as a function of J for different 
values of W. The results corresponds to the case x=10%, 
p=0.7%, and N=14. The Weiss mean field critical tempera- 
ture is plotted. 
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FIG. 5: Mn ions spin polarization as a function of T as ob- 
tained from Monte Carlo simulations and from Weiss mean 
field theory. The results correspond to x = 10%, p = 0.7%, 
J=W=3t. The MC results have been obtained in a cell of 
size N=14. 



the parameter W allows us to modulate the hybridization 
between the Mn ions orbitals and the host semiconduc- 
tor band. In the case of J—3t, we believe that a lower 
impurity band is not formed for any positive value of W. 

In Fig. [|we plot the critical temperature obtained from 
MC simulations as a function of J for different values of 
W, and for concentrations p=0.7% and x=10%. We have 
also plotted the critical temperature, Eq. 1^, obtained in 
the Weiss mean field theory @ The Weiss mean 
field T c increases quadratically with J and, as explained 
before, the approximation treats the magnetic impurities 
in the virtual crystal approximation, leading to a T^ CA 
independent of W. From the MC results, we observe 
that the MC critical temperatures are smaller than the 
T^ c s only for values of W very close to J. This is be- 
cause, for W=J, the effective impurity potential is zero, 
and the virtual crystal approximation for obtaining the 
Curie temperature seems to work well. In this case, W=J 
and small J, the thermal fluctuation correction to T^ CA 
is similar to that occurring in the standard Heisenberg 
model, ~ 25%. This becomes clear in Fig. |[ where we 
plot M(T) as obtained in MC simulation and compare it 
with the predictions of Weiss mean field. Both curves are 
very similar at low T, and only near the Curie tempera- 
ture do they separate, due to thermal fluctuation effects. 
For small values of W, the MC critical temperatures are 
larger than the TY CA, s. This is due to the charge located 
on the Mn ions which increases the effective magnetic 
field felt by the Mn spins. 

If the effective impurity potential J-W increases, holes 
would eventually form a localized impurity band, leading 
to a decrease of T c . As commented above, we do not get 
close enough to this limit, whose MC description would 
require a large basis beyond our present computational 
capabilities. 

We have also studied the dependence of the critical 
temperature on the Mn concentration, x. This is shown 
in Fig. ^|, where we plot T c as a function of x for J=3t, 
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FIG. 6: Curie temperature as a function of the Mn ion's con- 
centration as obtained from Monte Carlo simulations. The 
results correspond to p = 0.7%, J=3t and W=2t. The MC 
simulations have been done in a cell of size 7V=14. 



W=2t, and p=0.7%. Within the error margin, the criti- 
cal temperature increases linearly with expected for 
low impurity concentration. 

Finally, we have studied the dependence of T c on the 
hole charge density p. In Fig. we plot the MC results 
corresponding to the case J=3t, x — 10%, and two differ- 
ent values of W . For comparison, we have also plotted the 
Weiss mean field result. The Weiss mean field predicts 
a p 1 / 3 dependence, at low hole density. The MC results 
could be fitted perfectly to this dependence. Again we 
obtain that the Weiss mean field T c 's are higher that the 
MC results only for values of W near the value of the 
exchange coupling J. 

Before finishing our discussion on the MC results, a 
comparison with previous simulations is in order. Re- 
cently, Schliemann and coworkers [Esl have presented 
some MC results for (III,Mn)V semiconductors. They 
work in a model similar to ours, but with the site energy 
difference set to zero, W=0. They describe the motion of 
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FIG. 7: Curie temperature as a function of the hole concentra- 
tion, p, as obtained from Monte Carlo simulations and from 
Weiss mean field theory. The results correspond to x — 10% 
and J=3t. The MC simulations have been done in cells of 
size N=12, N=U, and N=16. 



FIG. 8: CPA critical temperature as a function of hole con- 
centration for x = 0.1, J = 3t, and different values of W . The 
Weiss (VCA) mean field critical temperature is also plotted 
for comparison. 



the carriers using a plane wave basis with a cutoff k c , and 
the exchange potential created by the magnetic impuri- 
ties is modeled with a Gaussian of thickness ao- They ob- 
tain that the critical temperature is substantially smaller 
than the estimated using Weiss mean field theory. This 
is opposite to our results: for W=0 our MC T c is higher 
than the mean field value. As the exchange constant and 
the densities (Mn and carriers) used in Ref. [Q are sim- 
ilar to ours, there is a clear contradiction between our 
results and those of reference p8| . A possible explana- 
tion for this discrepancy can be found in the wavelength 
cutoff used in rcf(^|), considerable larger than the im- 
purity Gaussian width, and, perhaps, unable to trap car- 
riers charge. If this were the case, their results could best 
correspond to our calculation for W—J, for which T c is 
smaller than TY CA . In reference the wavelength 

cutoff is near 25A, whereas the value of ao is smaller than 
5A. 

A comparison with experimental results is strongly de- 
pendent on the values of the parameters. For instance, 
in the case of t=0.85eV, J=M, x=W%, and p=0.7%, the 
critical temperature varies between 80K for W—J, and 
350K for Vy=0. This range would agree rather well with 
experimental data, assuming values of W near J. Given 
the oversimplified nature of our model, this comparison 
with experiments cannot be taken too far away. In any 
case, our results indicate clearly the relevant role of value 
of W on the magnetic transition: decreasing its value, if 
possible experimentally, would increase considerably the 
Curie temperature. 



B. CPA results. 

In this section, we present results for the phase dia- 
gram obtained with the CPA method explained in IV 



CPA results are parameterized by the density of states 
of the bare system, chosen to be that of the simple cu- 
bic lattice. Although unnecessary from the CPA point of 
view, we have assumed degenerate fermions. This is both 
physically correct for the relevant range of parameters, 
and makes the comparison with MC results more mean- 
ingful. We have verified that the inclusion of fermion 
temperature does not change results in any significant 
way. 

In fig. ^| we present the critical temperature for J = it 
as a function of carrier concentration in the relevant re- 
gion, for several values of W . As in the MC calculations 
(see Fig. |]), the most important effect of W is the de- 
pression of T c . Notice, though, that the overall temper- 
ature scale of the CPA calculation is lower than that of 
the MC results (compare with Fig. Q). This reverses the 
trend found in most statistical mechanics models, where 
the MC (exact) transition temperatures are higher that 
the mean field counterparts. We ignore whether this ef- 
fect is a genuine peculiarity of our model or rather the 
outcome of approximations. We observe that T c follows 
the p 1 / 3 dependence with carrier concentration expected 
to hold close to a 3-d band in the Weiss mean field ap- 
proximation, but this agreement between CPA and VCA 
is restricted to the small carrier concentration region (see 
below) . In addition, the very apparent effect of W on T c 
is completely missed in the VCA, as mentioned before. 

T c as a function of Hund coupling J is presented in 
fig. for several values of W . The prominent effect 
of W in reducing the critical temperature is very clear. 
Comparing with the MC data of Fig. |^, we see qualita- 
tively similar curves, but, again, the MC temperatures 
are higher than the CPA ones. It is worth comparing the 
CPA results of fig. | with those of the VCA (Eq. [jj]), 
which predicts a dependence of the form: TY° A oc xJ 2 , 
irrespective of the size of J and x. In the CPA treat- 
ment, this dependency only holds for small values of both 
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FIG. 9: CPA critical temperature as a function of J for 
x — 0.1 and p — 0.007. The different curves correspond to 
the following values of W (top to bottom): W = Ot, It, 2t, 
and 3t. Inset: Critical temperature as a function of impurity 
concentration for J = 3t, p = 0.007, and the same sequence 
of values of W (top to bottom) . 
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FIG. 10: CPA critical temperature as a function of hole 
density for J — 3t and the following values of W: W = Ot 
(continuous line), W = It (dashed line), and W — It (dashed- 
dotted line). Insets: Enlarged view of the phase diagrams 
showing the regions of phase separation for J = 3t, with W = 
Ot (left inset) and W = 2t (right inset). 



J and x, and significant deviations from the J 2 law are 
already evident in fig. ^, particularly for W = 0. We 
have also studied the dependence of T c with impurity 
density, as shown in the inset of fig. |^ for small values of 
x. As in the MC calculations, an approximately linear 
regime is obtained in this limit. This linearity suggests 
that the effective magnetic interaction between core-spins 
could be well approximated by a two-body interaction in 
this regime. Notice that, in spite of the VCA prediction 
(Eq. [jj]), this linear dependence does not hold for higher 
impurity density, or even in the concentration range of 
fig. ^| for W — 0, where non linear effects are already 
explicit. This behavior with J and x has been observed 
before for the present problem in a CPA calculation with 
W = and a model semicircular density of states p5[ . 
Its survival upon inclusion of disorder W ^ and for 
the real density of states (sc lattice) confirms that it is a 
robust feature of the model. 

It has been argued earlier that the effect of W on T c 
can be correlated with changes in carrier density on the 
impurity sites. Being the magnetic coupling mediated by 
mobile carriers, decreasing their density on the impurity 
sites (with the repulsive term W) should lead to a de- 
crease in T c . We can test this behavior analyzing the 
limit of an almost filled band, with electrons as carriers 
(remember that our bare fermions are holes). The ef- 
fect of W would then be the opposite and, therefore, T c 
should increase. To verify this idea while gaining a better 
perspective of the problem, we have plotted in fig. [l(] the 
phase diagram in the whole range of carrier concentra- 
tion, for three values of W . For W = 0, the problem (T c ) 
is particle-hole symmetric, and the ferromagnetic insta- 
bility is restricted to the band edges. Although quanti- 
tatively far from the small J limit, this confinement to 
band edges is qualitatively similar to the RKKY regime, 
where the same behavior can be traced back to the oscil- 



latory nature of the magnetic coupling. As anticipated, 
fig. |l0| confirms the asymmetric effect of W on both types 
of carriers: it decreases T c at the lower edge (holes) and 
increases T c at the upper edge (electrons). Notice that, 
for electrons, W acts as an attractive impurity, increas- 
ing the carrier density in the impurity sites, therefore 
increasing T c . It is worth mentioning that fig. ^ also 
contains the CPA results for W < 0: hole density p and 
W > is equivalent to having chosen W < and hole 
density 2 — p (i.e. electron density p). From this point of 
view, the physics close to the upper end of fig. |l^ upon in- 
creasing the value of W is controlled by the development 
of a magnetic impurity band, with magnetic behavior 
becoming closer to that of the Double-Exchange model 
|45f . Another interesting effect of increasing W shows 
up in fig. [l^: the spread of the magnetic instability from 
the lower edge into the band. We cannot offer a simple 
explanation for this fact. 

Finally, we have investigated the phase diagram look- 
ing for the presence of phase separation (PS). PS is a 
common feature in magnetic models with itinerant car- 
riers |5(], |3l], It seems to appear close to 
the magnetic transition temperature, when T c changes 
abruptly with carrier concentration p5| . In order to 
search for PS, we explore the existence of global minima 
in the thermodynamic potential, in addition to the usual 
local minimum employed to obtain the curve of T c . We 
observe small regions, close to the band edges, where the 
systems is unstable versus PS. This is shown in the in- 
sets of fig. [n] for two values of W. Comparing results for 
W = and W = 2t, we see that the PS region decreases 
with increasing W, due to the associated reduction in 
the dependence of T c with carrier concentration. In the 
limit of small electron concentration, similar PS areas 
(not shown here) are found with opposite behavior with 
W. 

As a final remark, we note that our calculation be- 
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comes qualitatively similar to a recent CPA analysis [E5l , 
carried out in the absence of impurity potential: W = 0. 
Quantitatively, however, it is different (even for W — 0) 
because of our using of the real density of states (sc lat- 
tice) as opposed to a model semicircular one. Never- 
theless, close to the band edges, the differences between 
both calculations for W = amount to a different choice 
of lattice regularization, or, equivalently, a rescaling of 
the Hund coupling J. 

VI. SUMMARY. 

We have studied a simple model for diluted magnetic 
semiconductors. In this model the carriers motion is de- 
scribed by a tight binding Hamiltonian and the hole-Mn 
ion coupling has two contributions: an antiferromagnetic 
exchange interaction between the Mn spin and the car- 
rier spin, J, and an electrostatic interaction, W , between 
the carriers charge and the potential arising from the Mn 
dopants. We have studied this model performing Monte 
Carlo simulations and mean field CPA calculations. Our 
main results are the following: 

1) The hole charge on the Mn ions has a relevant effect on 



the Curie temperature. The effective impurity potential 
J — W controls the hybridization of the Mn impurities 
with the host bands and the hole charge on the impu- 
rities. By increasing J — W , if possible experimentally, 
the charge on the Mn ions and the critical temperature 
would increase considerably. 

2) The Monte Carlo critical temperature is always higher 
than the corresponding CPA temperature. We also find 
that the Monte Carlo critical temperature is in general 
higher than the obtained from the simplest Weiss mean 
field theory. Only for values of J — W near cero the Weiss 
T c is higher than the Monte Carlo one. 

3) The CPA calculation shows the existence of small 
pockets of phase segregation instability close to band 
edges. 
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